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P.  L.  Rothwell,1  M.  B.  Silevitch,2  L.  P.  Block,3  and  P.  Tanskanen4 

A  model  of  ihe  westward  traveling  surge  (WTS)  and  the  generation  of  Pi  2  pulsations  is  presented  here. 
Previous  work  concentrated  on  the  motion  of  the  WTS  as  a  function  of  the  precipitating  electron  energy 
and  the  concurrent  generation  of  Pi  2  pulsations  via  a  feedback  instability.  Now  we  look  in  more  detail 
at  the  physical  assumptions  used  in  deriving  the  present  model  and  the  relations  between  the  zero-order 
and  the  first-order  solutions.  Constraints  are  placed  on  the  electron  temperature  asymmetry  in  the 
plasma  sheet  by  requiring  the  Pi  2  pulsations  to  be  bounded.  It  is  found  that  the  electron  temperature 
anisotropy  in  the  plasma  sheet  plays  a  major  role  in  determining  the  direction  in  which  the  surge  will 
propagate.  Narrower  surges  require  greater  electron  heating  parallel  to  the  magnetic  field  for  poleward 
motion.  More  energetic  electron  precipitation  is  predicted  to  produce  higher-frequency  Pi  2  pulsations. 
Pulsations  occur  in  multiple  bursts  with  the  time  interval  between  bursts  being  shorter  for  shorter  field 
lines.  Initial  amplitude  and  phase  conditions  are  cruc.al  in  determining  the  pulse  shape.  The  dominant 
period  of  the  Pt  2  pulsation  is  found  to  be  equal  to  twice  the  north-south  dimension  of  the  surge  divided 
by  a  term  which  is  proportional  to  the  poleward  velocity  of  the  boundary.  Finally,  we  show  that  the 
poleward  surge  velocities  and  Pi  2  pulsation  periods  as  measured  during  the  magnetospheric  substorm  of 
June  23,  1979.  are  consistent  with  our  model.  By  noting  the  direction  of  the  surge  motion,  one  can  use 
the  model  to  estimate  the  magnitude  of  the  polarization  electric  field.  We  find  that  it  is  consistent  with 
zero  for  the  onsets  considered. 
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1.  Introduction 

The  westward  traveling  surge  (WTS)  is  a  large  region  of 
auroral  brightening  that  occurs  near  local  midnight  during 
substorm  onsets.  This  region  generally  moves  in  a  northwest¬ 
ward  direction  but  at  times  is  seen  to  move  even  eastward. 
The  WTS  has  also  been  identified  as  the  source  of  Pi  2  pulsa¬ 
tions.  What  we  have  tried  to  do  over  the  past  few  years  is  to 
develop  a  unified  model  for  the  WTS  that  explains  both  the 
motion  of  the  surge  and  the  generation  of  the  Pi  2  pulsations 
[Rothwell  et  ul..  1984,  1986]  (hereinafter  referred  to  as  paper  1 
and  paper  2.  respectively).  Our  approach  has  been  to  model 
the  ionospheric  response  to  the  precipitating  electrons,  rather 
than  first  identifying  a  promising  magnetospheric  mechanism 
and  then  determining  its  effect  upon  the  ionosphere.  While 
this  may  appear  to  be  a  somewhat  backward  approach,  it  has 
the  advantage  of  starting  with  a  tractable  portion  of  the  prob¬ 
lem  whose  solution  imposes  conditions  on  the  far  less  trac¬ 
table  magnetospheric  source. 

It  is  implicitly  assumed  in  our  model  that  the  Inhester- 
Baumjohann  current  model  [Inhesler  et  a!..  1981;  Baumjo- 
hann.  1983]  represents  the  ionospheric  currents  inside  the 
surge  region  (see  Figure  1).  It  is  the  motion  of  this  current 
system  and  its  coupling  to  the  magnetosphere  that  are  the 
subject  of  this  paper.  While  the  recent  Viking  results  indicate 
that  the  formation  of  a  surge  may  be  far  more  complex  than 
previously  thought  [ Rostoker  et  al„  1987],  the  present  model 
is  still  considered  applicable  for  individual  "hot  spots"  as  well 
as  for  the  more  classical  type  of  surge  formation.  Note,  how- 
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ever,  that  our  model  deals  with  the  ionospheric  response  to 
electron  precipitation,  so  that  direct  changes  in  the  mag-  l 
netosphcric  sourcets).  such  as  the  creation  of  a  new  hot  spot  or 
the  fading  of  an  old  one,  could  produce  dynamic  effects  not 
explainable  by  the  present  model. 

The  surge  is  created  by  an  external  electric  field  E0  (Figure  i 
I).  This  field  drives  a  westward  Pedersen  current  and  a  pole- 
ward  Hall  current  The  Hall  current  is  closed  off  into  the 
magnetosphere  by  precipitating  electrons  along  the  poleward 
WTS  boundary.  The  parameter  x  is  a  measure  of  current  clo-  | 
sure  efficiency,  Lysak  [1986]  has  discussed  the  physical  inter¬ 
pretation  of  the  x  parameter  in  terms  of  a  possible  mag-  , 
netospheric  source  and  its  coupling  to  the  flux  tube 
ionosphere  system.  The  reader  is  referred  to  that  paper  for 
more  details.  In  our  context,  full  closure  (x  =  1 )  implies  the  full 
continuation  of  the  ionospheric  Hall  current  into  the  mag¬ 
netosphere  via  field-aligned  currents.  In  paper  1  we  found  that 
the  WTS  motion  in  the  midnight  sector  is  controlled  by  (1)  the 
energy  and  flux  of  the  precipitating  electrons.  (2)  the  electron- 
ion  recombination  rate,  and  (3)  the  degree  of  current  closure 
on  the  poleward  boundary  of  the  surge. 

On  the  basis  of  work  by  Rostoker  and  Samson  [1981], 
Samson  [1982].  Samson  and  Rostoker  [1983],  Pashin  el  al. 
[1982],  Lester  el  al.  [1984],  Singer  et  al.  [1983,  1985],  and 
Gelpi  et  al.  [1987],  there  is  clearly  an  observed  relationship 
between  the  WTS  and  the  generation  of  Pi  2  pulsations  during 
substorm  onsets.  It  is  natural  therefore  to  look  for  ways  in 
which  this  relationship  could  occur  using  the  current  system 
shown  in  Figure  I.  It  was  determined  in  paper  2  that  pertur¬ 
bations  in  the  north-south  current  in  the  WTS  could  produce 
standing  waves  due  to  the  reflection  from  conductivity  gradi¬ 
ents  along  the  surge  boundaries.  The  surge  serves  as  an  ac 
port  to  an  equivalent  transmission  line  [Soto,  1982]  formed  by 
the  attached  magnetic  field  lines.  Under  the  right  conditions 
this  port  can  resonate  with  the  Alfven  waves  that  propagate 
along  the  equivalent  transmission  line,  and  Pi  2  pulsations  will 
be  generated  from  the  ionosphere  into  the  magnetosphere. 

One  basic  idea  in  our  model  is  that  in  addition  to  the  ac 
(Alfven  wave)  component,  there  is  a  quasi  steady  state  (dc) 
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Fig.  I  Inhester-Baumjohann  ionospheric  current  system  [/n- 
hesier  et  til..  1981;  Baumjohanr.,  1983]  for  the  westward  traveling 
surge  This  model  forms  the  basis  for  the  model  presented  in  this 
paper 
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component  of  electron  precipitation  that  is  dependent  upon 
the  temperature  and  density  of  the  plasma  sheet  as  well  as  on 
the  potential  difference  along  the  field  line.  This  dc  component 
modulates  the  ionospheric  conductivity  and  controls  the  speed 
and  direction  of  the  surge  motion  (paper  1). 

The  interior  of  the  surge  has  a  higher  conductivity  than  the 
exterior  because  of  enhanced  precipitation.  The  current  car¬ 
ried  by  precipitating  energetic  electrons  in  the  interior  surge 
region  is  assumed  to  be  precisely  balanced  by  upward  flowing 
lower-energy  ionospheric  electrons.  Therefore  ionospheric  cur¬ 
rent  closure  into  the  magnetosphere  is  assumed  to  occur  only 
at  the  surge  boundaries. 

During  substorm  onsets  there  is  a  transient  injection  of  elec¬ 
trons  into  the  ionosphere  from  the  plasma  sheet  The  associ¬ 
ated  transient  current  is  assumed  in  our  model  to  coincide 
with  a  transverse  Alfven  wave  [Baumjohann  and  Glassmeier. 
1984]  that  initiates  a  feedback  instability  between  the  mag¬ 
netosphere  and  the  ionosphere  in  the  surge  region.  This  insta¬ 
bility  produces  a  complicated  first-order  ac  component  in  the 
precipitation  current  that  rides  on  top  of  the  dc  component  as 
discussed  above.  The  period  of  these  oscillations  tends  to  fall 
within  the  Pi  2  band  (40-160  s).  Now  the  first-  and  zero-order 
effects  are  not  decoupled.  The  dc  precipitation  flux  and  energy 
control  the  frequency  and  damping  rate  of  the  Pi  2  pulsations, 
as  we  will  see  below.  This  approach  is  significantly  different 
from  that  of  Lysak  [1986],  which  assumes  that  all  the  electron 
*  precipitation  is  associated  with  the  Alfven  wave.  Also,  our 
model  considers  the  ionospheric  response  to  the  Alfven  wave 
to  be  localized  in  the  surge  region,  not  global  in  nature,  as  in 
‘  the  work  of  Zhu  and  Kan  [1987],  Comparisons  of  our  work 
with  that  of  Kan  et  al.  [1984],  Kan  and  Kamide  [1985],  and 
Kan  and  Sun  [1985]  are  more  fully  discussed  in  section  6. 

In  the  work  of  Rothwell  et  al.  [1988]  the  results  of  Fridman 
and  Lemaire  [1980]  were  used  to  relate  the  dc  precipitation 
flux  and  energy  to  the  electron  temperature  anisotropy  in  the 
plasma  sheet  and  to  the  potential  drop  along  the  field  line.  It 
was  found  that  the  damping  rate  of  the  Pi  2  pulsations  as 
calculated  from  the  results  of  Fridman  and  Lemaire  closely 
followed  the  envelope  of  the  grow  th  curves  as  calculated  from 
the  feedback  instability  theory  of  Sato  [1982],  The  physical 
constraint  that  the  Pi  2  pulsations  be  damped  imposes  a  mini¬ 
mum  allowable  value  on  the  ratio  /,/a,  where  /,  is  the  north- 
south  dimension  of  the  surge  and  a  is  the  closure  parameter, 
as  defined  above.  It  was  found  that  narrower  surge  regions 
( —  100  km)  required  preferential  electron  heating  parallel  to  B 
(the  magnetic  field)  in  order  for  poleward  motion  to  take 
place.  For  thicker  surge  regions  ( —  300  km),  poleward  motion 
was  allowed  even  when  the  electron  temperature  perpendicu¬ 
lar  to  B  dominated  the  parallel  temperature.  The  precise  value 
of  the  threshold  ratio,  which  we  denote  by  (/,/«)„  depends  on 
the  electron  temperature  ratio  Tj/7j|  in  the  plasma  sheet.  In 
this  way  the  dynamics  of  the  surge  motion  is  controlled  by  the 


electron  temperature  anisotropy  in  the  plasma  sheet  bound¬ 
ary. 

In  section  2  we  will  review  the  zero-order  WTS  theory  as 
presented  in  paper  1.  In  section  3.  Pi  2  pulsations  will  be 
treated  according  to  the  first-order  theory  developed  in  paper 
2.  Section  4  deals  with  the  pulse  shapes  and  their  dependence 
on  initial  conditions  and  on  the  electron  precipitation  energy 
Section  5  compares  the  model  with  three  onsets  observed  on 
June  23.  1979.  Finally,  section  6  summarizes  and  presents  the 
conclusions  of  the  paper.  In  the  appendix  we  have  constructed 
a  very  simple  convection  model  to  justify  the  assumptions 
used  in  interpreting  the  experimental  data. 

2.  Zero-Order  Equations 

In  this  section  we  briefly  review  the  motion  of  the  poleward 
boundary  as  derived  in  paper  I.  Along  the  boundary  is  a 
conductivity  gradient  that  must  be  consistent  with  the  electron 
precipitation  from  the  magnetosphere.  The  ionospheric  ioniza¬ 
tion  density  in  the  gradient  region  is  governed  by  the  conti¬ 
nuity  equation  as  given  by 

<\V  <'/  =  Qj{l  e  -  of.\2  III 

where  Q  is  the  energy-dependent  ionization  efficiency  as  given 
by  Rees  [1963],  ju  e  is  the  electron  precipitation  flux.  <i,  is  the 
electron-ion  recombination  rate,  and  N  is  the  ionization  den¬ 
sity  in  the  ionosphere.  The  current  closure  on  the  poleward 
boundary  is  given  by 

,/j(  =  -x  i'J  e\  (2) 

where  J  is  the  poleward  Hall  current  and  x  is  the  closure 
parameter  along  the  poleward  boundary.  We  use  a  coordinate 
system  in  which  v  points  north,  y  points  west,  and  r  points  to 
the  zenith  J  is  related  to  ,V  by  ./  =  £HE0.  where  LH  is  the 
height-integrated  Hall  conductivity.  Now  ~  e,\  B.  where  B 
is  the  magnetic  field  at  the  ionosphere. 

By  combining  these  relationships  with  (2).  equation  (I)  be 
comes  a  wave  equation  if  we  momentarily  ignore  the  recombi¬ 
nation  term.  The  phase  velocity  of  this  ionization  wave  is 
given  by 

lv  -  Qhl^  (3<il 

where  h  is  the  ionospheric  height  over  which  Q  is  significant 
and  lj  -  F.„  B  ^  0.25  km  s. 

A  solution  for  the  velocity  of  the  poleward  boundary  in¬ 
cluding  electron-ion  recombination  effects  has  also  been  ob¬ 
tained  in  paper  I.  Linder  the  assumption  of  constant  closure  it 
was  found  that  the  equations  greatly  simplified  by  transform¬ 
ing  the  system  to  a  coordinate  frame  moving  at  the  boundary 
speed.  (A  more  exact  treatment  of  boundary  propagation, 
however,  must  take  into  account  time-  and  space-dependent 
closure.)  In  the  moving  frame  the  conductivity  profile  of  the 
boundary  remained  constant  and  matched  the  Hall  conduc¬ 
tivity  inside  the  surge  at  \  =  0.  (Note  that  v  is  the  poleward 
coordinate  here,  rather  than  as  in  paper  1.)  In  the  stationary 
frame  the  poleward  boundary  speed  as  given  by  equation  (22) 
in  paper  I  is 

V  -  F,  -  <752, f0 2 £„*  :/j|ft  G  =  orB  eh  (36) 

where  is  the  upward  closure  current  just  where  the  bound¬ 
ary  joins  the  surge  interior  and  £„0  is  the  height-integrated 
Hall  conductivity  in  the  surge  interior,  which  is  maintained  by 
a  precipitating  flux  jw;e.  Using  £wo  =  ehN0,‘B.  where  N0  is 
the  zero-order  ionization  level  inside  the  surge  region  which  is 
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given  by  [Qjlu  ea,]1  we  have 

V  -  Lr[  I  -  R,R„]  (3c) 

where  /?,  =  QJ\\,  e  and  Rh  =  QkjHb.e  are  the  ionization  rates  in 
the  surge  interior  and  on  the  poleward  boundary,  respectively, 
and  are  approximately  proportional  to  the  energy  flux.  (Note 
that  Q:  and  Q„  are  the  respective  ionization  efficiencies,  which 
are  energy  dependent  [Rees,  1963].)  If  the  ionization  rate 
along  the  poleward  surge  boundary  is  greater  than  the  ioniza¬ 
tion  rate  that  is  required  to  sustain  the  zero-order  level  of 
ionization  (IV0)  in  the  surge  interior,  then  the  surge  moves 
poleward.  If  it  is  less,  the  surge  moves  equatorward,  and  the 
surge  is  stationary  if  the  two  ionization  rates  are  equal.  The 
model  therefore  predicts  an  auroral  brightening  along  the 
poleward  surge  boundary  that  is  coincident  with  the  poleward 
leaps. 

Note  that  in  this  treatment  the  conductivity  profile  has  been 
assumed  to  remain  stationary  in  the  moving  frame,  which,  in 
general,  will  not  happen.  A  nonstationary  conductivity  profile 
leads  to  a  time  dependence  in  V  which  more  accurately  re¬ 
flects  the  surge  motion  [ Sileritch  el  ai.  1984],  The  propaga¬ 
tion  speed  of  any  profile,  however,  should  increase  in  propor¬ 
tion  to  the  excess  of  the  ionization  rate  along  the  boundary 
over  that  of  the  surge  interior.  After  all.  that  is  what  controls 
the  propagation  of  the  ionization  wave.  Therefore  for  short 
periods  right  after  onset  the  measured  velocities  should  be  well 
represented  by  (3a)  and  (3c).  This  assumption  is  utilized  below 
in  comparing  the  model  with  data  for  the  three  substorm 
onsets  of  June  23.  1979,  when  data  for  the  time  evolution  of 
the  conductivity  gradient  are  unavailable 

3.  Firs i -Order  Equations 

It  is  well  documented  that  Pi  2  pulsations  are  fundamen¬ 
tally  related  to  auroral  breakup  and  substorm  onsets  [Suito. 
1961 :  Rosloker  and  Samson.  1981 :  Samson  and  Rostoker.  1983: 
Samson.  1982].  As  pointed  out  by  Samson  [1982],  the  Pi  2 
pulsations  occur  simultaneously  with  or  before  all  other  iono¬ 
spheric  phenomena  associated  with  breakup.  These  results  are 
consistent  with  those  of  Sint/er  el  ai  [1983.  1985]  and  Gelpi  el 
ai  [1987],  who  used  magnetometer  data  from  the  Air  Force 
Geophysics  Laboratory  (AFGI.)  magnetometer  chain.  Their 
results  strongly  imply  that  the  Pi  2  source  is  located  approxi¬ 
mately  I  hour  to  the  east  of  the  western  surge  edge.  Pashin  el 
ai  [1982]  carried  out  a  study  on  Pi  2  pulsations  during  the 
passage  of  the  WTS  under  three  successive  auroral  breakups. 
They  found  that  the  largest  Pi  2  pulsation  amplitudes  were 
colocated  in  the  region  of  the  brightest  auroras.  In  a  similar 
study,  Baranskiy  el  ai  [1980]  observed  that  the  highest  Pi  2 
frequencies  occurred  near  local  midnight.  Stuart  el  ai  [1977] 
noted  a  correlation  between  the  modulation  in  precipitating 
electron  flux  in  the  auroral  zone  and  the  coincident  Pi  2  pulsa¬ 
tions.  Maltsev  el  ai  [1974]  suggested  that  Pi  2  pulsations  were 
a  result  ot  the  brightening  of  the  aurora  that  led  to  an  injec¬ 
tion  of  an  Alfven  wave  from  the  ionosphere  into  the  mag¬ 
netosphere.  This  leads  to  field  line  oscillations  and  resonances 
in  the  Pi  2  frequency  range. 

On  the  basis  of  the  large  volume  of  evidence  that  Pi  2 
pulsations  and  substorm  onsets  are  related,  it  is  therefore  logi¬ 
cal  to  explore  whether  the  surge  model  developed  in  paper  1 
can  explain  Pi  2  pulsations.  We  do  this  by  looking  at  the 
first -order  expansion  of  the  basic  equations.  Equations  (I)  and 
(2)  can  be  expanded  into  zero-  and  first-order  terms.  Along  the 
boundary  the  zero-order  terms  reproduce  the  results  in  section 


2,  while  they  lead  to  a  null  result  in  the  interior  where  N0  is 
assumed  constant.  The  first-order  terms,  on  the  other  hand, 
are  defined  in  the  surge  interior,  where  the  analysis  of  Sato 
[1982]  can  be  applied.  They  include  the  effect  of  Alfven  waves 
propagating  along  the  attached  field  lines.  Their  effect  at  the 
ionosphere  is  seen  as  an  equivalent  ac  impedance  [ Salo,  1982] 
which  is  given  by 

7.  =  iZ„  cot  (<’>,!.  L,  —  nit  I  (4a| 

where  Z0  —  VA  is  the  characteristic  impedance  of  the  equiv¬ 
alent  transmission  line,  to,  is  the  frequency  of  the  Alfven  wave, 
n  is  the  mode  number,  L  is  the  length  of  the  field  line  between 
the  ionosphere  and  the  equator,  and  V A  is  the  Alfven  speed 
( =  1 000  km/s)  in  the  magnetosphere.  Equation  (4a)  assumes 
perfect  reflection  of  the  Alfven  wave  at  the  equator  (i.e..  zero 
magnetospheric  conductivity),  in  contrast  to  the  finite  mag- 
netospheric  conductivity  model  developed  by  Lysak  [1986], 
The  first-order  precipitation  current  carried  by  the  Alfven 
wave  is  related  to  Z  through  Z/j(  =  V,  •  £,.  Implicit  in  this 
equation  is  the  transmission  line  assumption  and  the  propaga¬ 
tion 'reflection  of  Alfven  waves  within  the  flux  tube  connecting 
the  ionosphere  to  the  magnetosphere.  Now,  as  shown  in  paper 
2.  this  relation  when  combined  with  the  first-order  expansion 
of  (1)  and  (2)  leads  to  dispersion  relations  for  the  frequency 
and  growth  rates  of  the  individual  modes. 

Wc  will  briefly  summarize  the  derivation  as  given  in  paper 
2.  with  special  care  being  taken  regarding  whether  the  quan¬ 
tities  are  defined  on  the  surge  boundary  or  in  the  surge  in¬ 
terior.  Following  paper  2.  the  first-order  continuity  equation 
in  the  surge  interior  gives  us 

<'N(  el  -  —  Qhjtl,  eNoh  -  2 *N,  (4/>l 


x  =  °,H0h  h 

M,  is  the  height-integrated  perturbed  ionization  density  nor¬ 
malized  to  the  height-integrated  zero-order  density,  and  /,,,■  is 
the  perturbed  (first-order)  component  of  the  field-aligned  cur¬ 
rent.  which  is  given  by 

Jili”  —J„VNii\  —  ZZH0)  (4c) 


the  surge  interior  and  is  considered  to  be  constant  there.  On 
the  boundary,  however,  J„  closes  off  into  the  magnetosphere, 
giving  rise  to  the  zero-order  solutions  derived  above.  Note 
that  a  is  the  closure  parameter  as  defined  on  the  surge  bound¬ 
ary  while  Zlw  is  defined  in  the  surge  interior.  The  first-order 
current  fluctuations  give  rise  to  no  net  current  flow  into  and. 
out  of  the  interior  region.  Inserting  (4c)  into  (46 1  leads  to 
defining  the  parameter  I',  =  %QihVA,  which  is  analogous  to' 
the  boundary  velocity  parameter  as  defined  in  (3a)  and  (3c) 
and  which  we  now  call  VTh.  Note  that  Ki  =  Ks  only  if  the 
energy  of  the  zero-order  precipitation  is  the  same  in  the  surge 
interior  as  it  is  on  the  boundary.  In  the  absence  of  measure¬ 
ments  of  the  electron  precipitation  energy  we  assume  that  the 
velocity  of  the  poleward  boundary  gives  a  reasonable  measure 
of  VfI.  This  assumption  is  used  below  in  section  5  and  dis¬ 
cussed  in  the  appendix. 

The  first-order  dispersion  relations  are  then  given  by 
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□ 
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Pt2  GROWTH  AND  DECAY  RATES  VS  V,/X, 

10  V./L- 0.005  A,/a-  130km  T,*04MV  T±  ■  0.2k«V 


Fig.  2.  Growth  and  damping  curves  for  Pi  2  pulsations  as  derived 
from  model  The  growth  curves  are  the  nearly  vertical  lines,  and  the 
damping  (decay)  curve  follows  the  envelope  of  the  growth  curves.  The 
tick  marks  on  the  decay  curve  map  the  precipitating  electron  energies 
to  the  abscissa.  The  abscissa  is  plotted  in  units  of  F  / where  i.x  is 
the  north-south  extent  of  the  surge.  In  this  example  we  took  /.x 
y  -  I  JO  km.  V,  L  —  0.005.  LpoZ„  -  10.  TJ,  -  0  4  keV.  and  T  -  0  2 
keV.  See  text  for  definitions. 

w  here  <  jr  and  <o,  are  the  real  and  imaginary  components  of  the 
frequency,  respectively.  V  is  equal  to  Z  as  defined  in  (4a)  times 
the  zero-order  Pedersen  conductivity  in  the  ionosphere. 

\,,k  is  the  height-integrated  equilibrium  ionization  density  as 
expected  from  the  zero-order  (dc)  electron  precipitation  inside 
the  surge.  From  now  on  we  set  N0lx  =  h\'n.  where  ,V0  is  the 
local  zero-order  ionization  density  at  the  maximum  ionization 
altitude  for  a  given  precipitation  energy.  (Note  that  in  equa¬ 
tion  1 18)  of  paper  2,  li  of  equation  (6)  above  is  missing  and  the 
definition  given  in  equation  (1 1)  of  paper  2  should  replace  x  in 
equations  (12)  and  (16)  of  paper  2.) 

Now  if  the  zero-order  electron  precipitation  energy  inside 
the  surge  is  of  the  order  of  a  couple  keV.  then  Vxi  is  approxi¬ 
mately  a  couple  kilometers  per  second  (paper  I).  In  addition,  if 
the  north-south  surge  dimensions  (z.x)  are  of  the  order  of 
hundreds  of  kilometers,  then  the  characteristic  frequency  (V,, 
2/.x)  for  ionospheric  waves  is  close  to  that  of  Pi  2  pulsations. 
In  order  for  a  resonance  to  occur  the  characteristic  iono¬ 
spheric  frequency  just  derived  must  be  greater  than  or  equal 
to  the  frequency  of  the  fundamental  toroidal  mode  for  the 
attached  field  lines.  Sini/er  et  al.  [1981]  found  that  the  fre¬ 
quency  of  the  fundamental  toroidal  mode  was  smallest  at  local 
midnight,  reaching  the  Pi  2  range  of  0.02  s'  1  at  A  =  66  and 
0.005  s  1  at  A  =  68  .  This  suggests  that  the  formation  of  a 
surge  at  auroral  latitudes  near  local  midnight  sets  up  a  reso¬ 
nance  condition  between  the  ionospheric  waves  and  the  mag¬ 
netic  field  line  oscillations  which  is  triggered  by  the  feedback 
instability,  as  described  here. 

There  are  both  reactive  (+30  and  capacitive  ( —  A")  solu¬ 
tions  to  (5)  and  (6).  The  reactive  solutions  lead  to  growing  Pi  2 
pulsations  from  the  feedback  instability,  while  the  capacitive 
solutions  quench  the  instability.  Reference  is  made  to  Sato 
[1982]  and  references  therein  regarding  the  details  of  these 
effects.  Samson  [1982]  found  that  the  Pi  2  oscillations  origi¬ 
nating  from  field-aligned  currents  are  phase-shifted  with  re¬ 
spect  to  Pi  2  pulsations  arising  from  the  electrojet.  In  the 
present  model  this  effect  comes  from  the  inductive  nature  of 
the  field  lines  which  allows  growing  solutions. 


Note  that  the  first-order  dispersion  relations  also  depend  on 
the  zero-order  precipitation  energy  (i.e..  the  field-aligned  po¬ 
tential  drop)  through  Vxl.  We  use  the  theory  of  Fridman  and 
Lemaire  [1980]  to  estimate  the  magnitude  of  the  zero-order 
precipitation  for  various  field-aligned  potential  differences  and 
plasma  sheet  parameters. 

Figure  2  shows  a  separate  plot  of  the  Pi  2  growth  and  decay 
terms  from  (6)  as  a  function  of  lx  zx.  The  growth  curves  are 
the  nearly  vertical  lines  that  tend  to  cross  over  each  other.  The 
damping  curve  is  the  nearly  straight  line  that  closely  follows 
the  envelope  of  the  growth  curves.  The  damping  curve  was 
calculated  according  to  the  Fridman  and  Lemaire  [1980] 
theory  with  electron  plasma  sheet  temperatures  of  TJ,  =  0.4 
keV  and  T  =  0.2  keV  and  with  the  density  of  the  plasma 
sheet  assumed  to  be  0.3  cm  “  J.  The  tick  marks  denote  the 
locations  of  various  precipitation  energies.  Note  that  more 
energetic  precipitation  turns  on  higher-frequency  modes.  For 
Figure  2  we  used  Lro  =  10  mhos.  Z0  =  1  ohm.  T,  L  =  0.005 
s'1,  and  zx  x  =  I  50  km. 

If  one  inserts  (5)  into  (6)  and.  for  the  moment,  ignores  the 
decay  term,  then  one  can  find  the  maximum  growth  rate  as  a 
function  of  A  .  It  i>  <  >r\  =  r>r  =  nl  \  /.x  [ Rotltwell  et  al..  1988] 
The  condition  that  the  Pi  2  pulsations  be  damped  (i.e  (hat 
2it,.V(,  2  itl>  zj  leads  to  a  minimum  allowable  threshold  on 

x.  w  hich  to  a  high  degree  of  accuracy  is  independent  of  the 
zero-order  precipitation  energy  This  threshold  depends  on  the 
plasma  sheet  parameters  through  the  Fridman  and  Lemaire 
[1980]  relation  We  again  set  the  plasma  sheet  number  density 
to  0.3  cm  1  and  calculated  the  threshold  value  of  /.x  x  Ky 
assuming  a  linear  fit  to  the  precipitation  current  for  field- 
aligned  potential  drops  between  5  kV  and  50  kV.  Figure  3 
shows  a  plot  of  the  threshold  values  (z.x  x),  versus  the  electron 
temperature  anisotropy  T  7],  in  (he  plasma  sheet.  An  analyti¬ 
cal  expression  for  the  dependence  of  the  threshold  values  on 
T  .  7J|.  and  (he  plasma  sheet  number  density  .V  can  be  found. 
This  is  done  by  noting  lhat  \„  >.  fl>.  which  leads  to  (z.x  x),  > 

I  1  :  ('/j,1  *Npx'  ’).  where  <I>  is  the  field-aligned  potential  drop 

There  is  a  simple  physical  explanation  for  associating  an 
enhanced  7],  with  poleward  motion  of  the  conductivity  slab 
An  increase  in  7J(  causes  an  increase  in  jnt,  [ Fridman  and  Le¬ 
maire.  1980],  This  increases  the  degree  of  modification  of  the 
conductivity  gradient  along  the  poleward  boundary,  which 


Fig.  5.  The  damping  threshold  for  Pi  2  pulsations  as  a  function  of 
the  electron  temperature  anisotropy  in  the  plasma  sheet  using  the 
model  of  Fridman  and  Lemaire  [1980]  for  the  precipitating  electron 
flux.  When  the  measured  value  of  /t  x  exceeds  the  threshold  value, 
the  Pi  2  pulsations  are  damped 
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causes  a  faster  poleward  surge  expansion.  The  net  effect  is  that 
anything  that  increases  j^b  along  the  poleward  boundary  will 
cause  a  more  northward  motion  of  the  surge  region.  An  en¬ 
hancement  of  juh  also  reduces  the  north-south  polarization 
field  (see  Figure  1).  which  decreases  the  westward  Hall  current. 
Now  an  enhancement  in  temperature  anisotropy  (if  it  is 
global)  will  tend  to  impose  an  increase  in  jn  at  the  surge  head. 
In  our  simplified  steady  state  model  this  must  be  offset  by  a 
compensating  decrease  in  the  field-aligned  potential  [ Fridman 
and  Lemaire,  1980]  at  the  surge  head  in  order  for  jn  there  to  be 
continuous  with  the  diminished  westward  ionospheric  current. 
Our  model  therefore  also  suggests  that  enhanced  7^  is  accom¬ 
panied  by  a  brightening  along  the  poleward  boundary  and  a 
dimming  in  the  region  near  the  surge  head. 

The  opposite  argument  holds  for  preferential  electron  heat¬ 
ing  perpendicular  to  the  magnetic  field  in  the  plasma  sheet.  In 
that  case  the  westward  Hall  current  due  to  the  polarization 
electric  field  is  increased,  and  therefore  the  field-aligned  poten¬ 
tial  at  the  surge  head  must  increase  to  maintain  current  conti¬ 
nuity  [ Fridman  and  Lemaire.  1980].  These  effects  favor  a  more 
westward  motion  of  the  surge  region.  The  corresponding  argu¬ 
ments  suggest  that  enhanced  TL  will  also  cause  a  dimming 
along  the  poleward  boundary  accompanied  by  a  brightening 
near  the  surge  head. 

In  summary,  according  to  our  model,  poleward  surge 
motion  is  associated  with  electron  heating  in  the  plasma  sheet 
parallel  to  the  magnetic  field.  This  heating  could  be  caused  by 
Fermi  acceleration  as  the  tail  field  becomes  more  dipolar. 
Preferential  perpendicular  heating  from  conservation  of  the 
first  invariant  is  associated  with  westward  surge  motion. 

What  happens  when  the  surge  size  is  below  the  threshold 
limit  with  x  =  I  ?  In  that  case  we  have  positive  growing  pulsa¬ 
tions.  The  period  of  the  wave  with  the  fastest  growth  rate,  and 
therefore  the  period  that  is  characteristic  of  the  Pi  2  pulsation 
is  2/.v  Ij.  In  the  limit  of  (3<t)  the  surge  is  expanding  at  a 
velocity  I'.,  which  means  the  surge  triples  in  size  in  one  pulsa¬ 
tion  period.  The  threshold  for  damping  is  quickly  reached,  and 
the  pulsation  has  a  finite  amplitude.  This  feature  validates  the 
linear  approximation  used  in  deriving  the  Pi  2  pulsations. 
However,  the  rapid  poleward  motions  mean  that  the  zero- 
order  parameter  is  changing  on  a  time  scale  comparable 
with  the  first-order  solutions.  At  the  very  least  this  leads  to  a 
time-dependent  frequency  which  is  outside  the  scope  of  our 
solutions.  This  problem  is  resolved  if  the  surge  moves  in  spurts 
or  jumps.  In  that  case  the  motion  of  the  boundary  induces  a 
standing  ionization  wave  inside  the  surge  with  a  phase  veloci¬ 
ty  which  we  assume  to  be  close  to  the  boundary  velocity. 
From  (5)  this  ionization  wave  is  in  resonance  with  the  fastest 
growing  Pi  2  mode  (X  =  I),  and  the  present  model  applies 
between  jumps.  The  assumption  is  that  the  jumps  occur  often 
enough  that  the  critical  surge  size  is  reached  in  a  timely 
manner  and  the  pulsation  amplitudes  are  not  unreasonably 
large.  A  parallel  resistance  along  the  field  line,  such  as  that 
measured  by  Weimer  et  al.  [1985.  1987],  would  also  contrib¬ 
ute  to  keeping  the  amplitudes  finite  (C.  K.  Goertz.  private 
communication,  1987).  During  jumps  a  nonlinear  treatment 
such  as  that  of  Lysak  [1986.  1985],  which  incorporates  the 
generation  of  PiB  pulsations,  is  probably  necessary.  The  ob¬ 
servations  of  Opgenoorth  et  al.  [1980]  and  Basinger  et  al. 
[1981]  suggest  that  PiB  type  pulsations  lasted  as  long  as  the 
local  onset-connected  field-aligned  currents  were  growing.  On 
the  other  hand,  if  the  multiplicative  factor  shown  in  (.V)  is 
sufficiently  small,  then  the  change  in  frequency  is  adiabatic, 
and  the  present  model  applies  in  its  entirety. 


fig.  4.  (Xvay  lime  of  the  individual  modes  shown  in  figure  2  for 
a  precipitating  eleelron  energy  of  5  keV.  The  dots  denote  specific 
modes  labeled  n.  Note  that  after  several  minutes  the  pulse  will  become 
quasi-monoehromalie  wnh  ihe  center  frequency  corresponding  lo  the 
fastest  growing  mode  This  is  consistent  with  ihe  dala  shown  in 
F  igure  9, 

From  Figure  3  we  see  that  smaller  damping  threshold 
values  for  /.x  result  when  the  plasma  sheet  electrons  are  prefer¬ 
entially  heated  parallel  to  the  magnetic  field.  Now  if  the  north- 
south  surge  dimension  is  below  the  threshold  value  for  a  given 
plasma  sheet  model,  then  x  must  be  less  than  I  in  order  for  the 
Pi  2  pulsations  to  be  damped.  In  paper  I  we  showed  that  the 
direction  of  the  surge  motion  is  highly  dependent  on  x  For  x 
values  less  than  I.  westward  motion  is  predicted,  while  x  g  1 
allows  poleward  or  even  eastward  motion.  Therefore  accord¬ 
ing  to  our  model  the  electron  temperature  anisotropy  in  the 
plasma  sheet  together  with  the  physical  requirement  that  the 
Pi  2  waves  be  damped  plays  an  important  role  in  determining 
which  way  the  current  system  of  Figure  I  is  going  to  move. 
The  result  is  that  embryonic  surge  regions  reutiire  substantial 
preferential  heating  parallel  to  B  for  polewaid  leaps  to  occur. 
Some  specific  pulse  shapes  will  now  be  exam.ned 

4.  Pi  2  Pit.se  Si  .res 

Recall  from  F  igure  2  that  the  net  decay  rate  of  each  mode  is 
the  difference  between  the  decay  curve  (2(Tr/V'„l  and  the  indi¬ 
vidual  growth  curves  at  a  specific  precipitation  energy.  Figure 
4  shows  the  calculated  decay  time  for  each  mode  for  an  inci¬ 
dent  energy  of  5  keV  ard  a  surge  size  of  113  km.  All  other 
parameter  values  are  a  in  Figure  2. 

Thirty  modes  arv  excited  in  this  particular  example,  and 
each  mode  is  repi  evented  by  a  dot.  Note  that  the  decay  time 
as  a  function  4  frequency  is  highly  peaked  at  the  frequency 
value  found  in  section  3  for  the  maximum  growth  rate.  What 
Figure  4  tells  us  is  that  the  higher  and  lower  frequencies 
quickly  become  damped  and  that  the  pulse  propagates  around 
a  narrower  frequency  hand  at  longer  times  The  period  of  this 
pulsation  at  the  peak  frequency  is  approximately  25  s.  which  is 
consistent  with  the  Pi  2  frequency  range. 
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I-'ijj.  5,  A  sample  Pi  2  pulse  shape  (ampliluile  in  arbitrary  unilsl  as 
a  function  oi  time  for  the  same  input  parameters  as  shown  in  Figure 
4  I  his  is  a  linear  super  position  of  all  the  modes  assuming  equal 
initial  amplitudes  and  phases. 

F"  re  5  is  a  time  history  of  this  pulse.  The  initial  con¬ 
ditions  for  this  particular  plot  were  chosen  such  that  all  modes 
were  assumed  to  he  excited  with  the  same  amplitude  and  at 
the  same  phase  The  predicted  time  history  as  seen  from 
t-igure  5  is  that  of  a  large,  initially  oxerdamped  pulse  followed 
by  series  of  "ringing"  pulses  of  about  .3  min  in  duration  and 
6  7  min  apart  The  time  interval  between  "rings"  becomes 
longer  as  the  ratio  T,  /.  becomes  smaller.  Assuming  T,  stays 
constant,  this  implies  that  period  between  rings  is  proportion¬ 
al  to  the  length  of  the  magnetic  field  line 

I  he  specific  shape  of  a  given  pulsation  therefore  depends 
very  strongly  on  the  initial  conditions  of  its  formation.  The 
present  model  does  not  specify  what  these  conditions  are.  but 
it  does  give  the  time  evolution  of  the  pulse  once  the  initial 
conditions  are  given 

Despite  this  uncertainty,  however,  one  can  speculate  on 
various  possibilities  We  redid  the  plot  in  Figure  5  using 
random  initial  phases  I  betw  een  0  and  360  I  for  each  of  the  30 
modes  Figure  6  shows  the  results.  Note  that  the  pulse  has  an 
extremely  complicated  shape  during  the  first  5  min  which 
completely  masks  the  underlying  sinusoidal  dependence.  In 
fact,  one  might  interpret  this  section  as  a  series  of  single 
pulses  The  randomness,  however,  destroys  the  precise  re¬ 
peatability  of  the  time  profile,  so  that  another  example  with 
the  same  parameter  values  would  look  somewhat  different.  It 
was  found  that  even  if  the  randomness  in  phase  was  as  little  as 
20  .  a  difference  with  Figure  6  was  not  discernible.  Given  the 
complexity  of  the  magnetosphere-ionosphere  system  and  the 
irregularity  of  the  driving  precipitating  flux,  some  incoherence 
in  mode  stimulation  is  to  be  expected.  This  incoherence  de¬ 
stroys  the  precise  cancellation  between  bursts,  as  seen  in 
Figure  5.  and  the  model  simulates  the  actual  data  more  close- 

13 

Randomization  of  the  initial  amplitudes  also  destroys  co¬ 
herence  and  produces  results  similar  to  those  seen  in  Figure  6 
for  random  initial  phases.  Therefore  Figure  6  is  representative 
of  a  wide  class  of  initial  conditions.  Now,  regardless  of  the 
initial  conditions,  the  selective  filtering  of  the  frequencies  as 
seen  in  Figure  4  leads  to  a  burstlike  type  of  behavior  after 
longer  times  Therefore  the  precise  nature  of  the  initial  con¬ 
ditions  seems  to  dominate  in  only  the  first  few  minutes  of  the 
pulse  train 


Figure  7  shows  a  coherent  pulse  train  for  an  incident  energy 
of  1.5  keV,  l'A  l.  —  0.001  s  '.  and  x  =  125  km.  Note  the 
same  general  features  as  in  Figure  5:  an  initial  overdamped 
pulse  followed  much  later  (30  min)  by  a  finite  pulse  train.  If 
this  plot  were  of  data,  it  would  take  a  very  astute  observer  to 
recognize  that  these  two  pulses  are  indeed  one  and  the  same. 
The  pulsation  burst  evolves  toward  the  mode  with  the  longest 
decay  time.  (See  the  above  discussion  in  regard  to  Figure  4  ) 
The  period  of  this  mode  is  given  by  Tp  =  2/,,/T',  and  in  this 
case  is  equal  to  122  s,  as  opposed  to  25  s  for  a  5-keV  electron 
precipitation  energy.  Therefore  one  testable  element  of  our 
model  is  that  the  frequency  of  the  later,  periodic  portions  of 
the  Pi  2  pulsations  should  increase  with  the  electron  precipi¬ 
tation  energy. 

Figure  8  shows  a  pulsation  with  the  same  input  values  as  in 
Figure  7  except  that  the  initial  phases  are  randomized  as  in 
Figure  6.  Note  that  the  horizontal  scale  is  approximately  one- 
half  that  of  F  igure  7.  Randomization  causes  the  various  modes 
not  to  cancel  and  produces  some  interesting  waveforms  during 
the  first  few  minutes.  In  fact,  the  first  two  peaks  resemble  a 
sequence  of  two  overdamped  single  pulses,  rather  than  a 
superposition  of  sinusoidal  modes.  The  underlying  sinusoidal 
structure  does  not  become  apparent  until  10  min  after  onset. 
The  pulsation  shape  is  close  to  the  measured  Pi  2  waveforms 
shown  in  paper  2  and  in  Figure  3  in  the  work  of  Southivnml 
ami  Smart  f  iyS0] 

It  should  be  pointed  out  that  the  precise  cancellation  of  the 
waveforms  between  bursts  seen  in  Figures  5  and  7  implies  a 
highly  systematic  structure  for  the  nonlinear  solutions  which 
were  obtained  numerically.  Investigation  is  proceeding  as  to 
what  this  structure  is  and  whether  it  can  be  expressed  in  an 
analytical  form  or  not. 

We  are  also  presently  looking  at  the  Fourier  transform  of 
these  time  profiles  in  order  to  make  comparisons  with  data. 
Owing  to  the  fact  that  the  decay  time  is  a  function  of  fre¬ 
quency  (see  Figure  4).  the  frequency  content  of  the  wave  is 
time  dependent.  This  means  that  the  Fourier  transform  is  sen¬ 
sitive  to  the  time  at  which  the  transform  is  taken.  Neverthe¬ 
less.  our  model  predicts  a  spiky  power  density  spectrum  with 
the  central  spike  located  at  the  frequency  corresponding  to 
minimum  damping. 


Fig.  6.  The  same  as  F  igure  5  but  now  the  linear  superposition  is 
carried  out  with  random  initial  phases.  Note  the  nonsinusoidal  be¬ 
havior  al  Ihe  beginning  This  figure  is  more  consistent  with  the  data 
than  Figure  5. 
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5.  Comparison  oe  Model  With 
Experimental  Data 

A  detailed  study  of  three  successive  substorm  onsets  on 
June  23,  1979,  was  carried  out  using  satellite-,  balloon-,  and 
ground-based  particle  and  field  data  [ Tanskanen  ei  at.,  1987], 
We  now  utilize  these  data  to  estimate  the  poleward  velocity  of 
the  surge  boundary  and  the  north-south  extension  of  the  surge 
for  each  of  the  three  onsets.  These  results  are  then  used  to 
estimate  the  period  of  the  fastest  growing  mode,  as  described 
above.  This  period  is  then  compared  with  the  measured  period 
of  Pi  2  pulsations  as  seen  by  ground-based  magnetometer 
stations. 

For  the  following  comparison  we  use  Table  3  of  Tanskanen 
et  al.  [1987],  The  first  onset  was  observed  at  Sodankyla  (lati¬ 
tude  67.4  .  longitude  26.6  ,  L  =  5.1)  at  2101  UT.  Thirty  sec¬ 
onds  later,  it  was  observed  at  Ivalo  (latitude  68.6  ,  longitude 
27.4  ,  l.  =  5.7).  and  30  s  beyond  that,  at  2102  UT,  it  was 
observed  at  Kevo  (latitude  69.8  ,  longitude  27.0  .  L  =  6.2). 
Now  Sodankyla,  Ivalo.  and  Kevo  are  close  together  on  the 
same  meridian,  so  that  the  successive  detection  of  the  onset 
gives  a  good  measure  of  the  initial  poleward  surge  velocity  V . 
We  find  V  -  4.4  km/'s  for  the  first  onset.  The  westward  com¬ 
ponent  of  the  surge  velocity  Vw  is  found  by  noting  that  the 
onset  occurs  at  the  balloon  S0622  (latitude  67.6  ,  longitude 
5.1  ,  L  ^  6)  10  min  after  it  reached  Sodankyla.  which  is  along 
the  same  geographic  latitude.  Therefore  Vw  =  1.6  km/s  for  the 
first  onset.  The  surge  is  estimated  to  include  Sodankyla  and 
Kevo.  which  gives  a  north-south  extension  ).x  of  about  280 
km.  For  a  complete  analysis  it  is  necessary  to  know  the  energy 
and  flux  of  the  electron  precipitation  both  inside  the  surge  and 
along  the  poleward  boundary.  Unfortunately,  this  level  of 
detail  is  not  available  However,  we  do  have  riometer  data 
which  give  an  indication  of  the  relative  intensity  of  the  precipi¬ 
tation  energy.  The  precipitation  was  most  intense  during  the 
first  onset.  Therefore  for  this  period  we  make  the  assumption 
that  the  ionization  rate  along  the  surge  boundary  exceeds  the 
ionization  rate  inside  the  surge  and  that  V  =  Fx(1  =  Fx,  =  4.4 
km  s.  (See  (3c).  the  discussion  following  (4c).  and  the  appendix 
for  more  details  regarding  these  assumptions )  In  this  case  the 
period  tTp  =  2/.x  Fxj)  of  the  fastest  growing  mode  is  127  s. 
Now  the  magnetometer  results  are  Fourier-analyzed  and  plot- 


Hg.  7  Time  history  of  a  Pi  2  pulsation  with  coherent  initial 
phases  and  1.5  keV  incident  energy.  Note  the  longer  period  in  com¬ 
parison  with  Figure  5.  The  higher  the  incident  energy  the  higher  the 
frequency  of  the  expected  pulsations. 


Fig.  8.  Time  history  of  a  Pi  2  pulsation  with  random  initial 
phases  and  1.5  keV  incident  energy.  This  pulse  profile  is  very  similar 
to  profiles  observed  in  the  data  (see,  for  example.  Rothwell  et  al. 
[1986]).  Note  once  again  that  the  pulsation  becomes  almost  mono¬ 
chromatic  after  about  10  min. 

ted  in  frequency  versus  time  contour  plots,  where  each  con¬ 
tour  represents  a  constant  amplitude.  The  contour  plot  of 
Kevo  is  shown  in  Figure  9,  with  the  onsets  denoted  by  short 
vertical  lines.  Just  after  the  first  two  onsets  the  maximum 
amplitude  (the  dense  band  of  concentric  ellipses)  occurred  at 
approximately  120  +  10  s,  with  a  secondary  peak  at  ~95  s.  A 
plot  of  the  Ivalo  magnetometer  data  shows  similar  results 
although  at  lower  amplitudes.  The  delay  of  the  monochro¬ 
maticity  relative  to  the  substorm  onsets  is  consistent  with  the 
model  results  shown  in  Figures  6  and  8.  For  the  first  onset  the 
delay  is  approximately  6-8  min.  Considering  the  assumptions 
involved  in  calculating  F„  and  the  excellent  agreement  (127 
s  versus  120  +  10  s)  may  be  somewhat  fortuitous.  However, 
the  overall  agreement  between  the  experimental  observations 
and  the  model  lends  a  great  deal  of  credence  to  the  model's 
basic  veracity. 

The  second  onset  of  Sodankyla  (2107:30  UT)  gives  similar 
results.  The  poleward  velocity  between  Sodankyla  and  Ivalo 
(2109  UT)  is  1.5  km/s.  and  between  Ivalo  and  Kevo  (2109:30 
UT)  it  is  4.4  km/s.  The  lower  velocity  between  Sodankyla  and 
Ivalo  is  explainable  by  the  fact  that  the  riometer  absorption  at 
Sodankyla  has  dropped  to  0.7  dB  from  1.3  dB  as  in  the  first 
onset.  This  implies  less  intense  precipitation  and  hence  a 
slower  speed  according  to  our  model  The  westward  velocity 
component  is  found  as  before  and  is  1.8  km/s.  We  take  the 
north-south  surge  dimensions  to  be  between  Sodankyla  and 
Kevo  as  before,  which  gives  286  km  in  geographic  coordi¬ 
nates.  Now  Andenes  is  geographically  south  but  magnetically 
north  of  Kevo  [see  Tanskanen  et  al.,  1987.  Table  3].  If  one 
prefers  the  surge  dimensions  in  magnetic  coordinates,  then  the 
surge  should  be  measured  between  Andenes  and  Sodankyla. 
In  that  case  we  find  a  distance  of  289  km.  The  period  of  the 
fastest  growing  mode  is  therefore  predicted  to  be  131  s.  Now 
from  Figure  9  we  see  that  the  period  with  the  maximum  am¬ 
plitude  is  around  127  s,  which  occurs  approximately  10  min 
after  the  second  onset.  The  agreement  is  again  excellent. 

The  third  onset  is  more  complicated  for  two  reasons.  The 
first  reason  is  that  the  poleward  boundary  now  passes  by  the 
stations  KA  0622  (latitude  70.4°,  longitude  —1.2  .  L  ~  8.5)  at 
2140  UT,  Bjornoya  (latitude  74.5  ,  longitude  19  ,  L  -  9.4)  at 
2143  UT.  and  HO  0623  (latitude  72.9",  longitude  7.6  .  L  ~  9.2) 
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Fiji  4.  A  contour  plot  of  the  littered  magnetometer  data  taken  at  Kevo  during  Ihe  period  of  three  substorm  onsets  on 
June  23.  1474  Each  contour  represents  a  given  amplitude  which  is  plotted  as  a  function  of  frequency.  The  abscissa  is 
universal  lime  in  units  of  minutes  starting  at  2058:00  UT  and  ending  at  2158:00  UT.  The  three  heavy  tick  marks  show  the 
time  of  the  three  onsets  The  more  concentrated  areas  of  concentric  rings  denote  regions  of  more  intense  wave  energy.  It  is 
shown  in  the  text  that  these  regions  are  numerically  very  consistent  with  the  expected  Pi  2  periods  as  derived  from  our 
model  Also,  the  lengthening  of  the  Pi  2  period  with  successive  onsets  is  consistent  with  the  poleward  expansion  of  the 
surge  as  predicted  by  the  model 


at  2145  UT.  These  stations  are  separated  in  both  latitude  and 
longitude,  which  complicates  the  determination  of  the  velocity 
components.  The  second  reason  is  that  the  third  onset  is  much 
less  intense  than  the  previous  two  onsets  according  to  the 
riometer  data.  This  implies  from  our  model  (see  (3c))  that  the 
measured  poleward  velocity  could  be  much  less  than  the  Ttl 
parameter  needed  to  estimate  the  associated  Pi  2  period.  We 
lind  that  the  measured  poleward  velocity  is  approximately 
equal  to  0.25  km  s  and  the  westward  velocity  is  2  3  kms. 
where  the  spread  is  due  to  whether  one  works  in  geographic 
or  magnetic  coordinates.  We  therefore  assume  for  the  third 
onset  that  the  precipitation  energy  is  the  same  inside  the  surge 
as  on  the  poleward  boundary  and  as  during  the  previous  two 
onsets.  According  to  this  scenario  the  surge  moves  because  of 
an  enhancement  of  the  boundary  precipitation  flux  of  a  few 
percent  over  the  precipitation  flux  inside  the  surge.  The  old 
value  of  4.4  km  s  for  U ,  will  therefore  be  used. 

The  north-south  extent  of  the  surge  during  the  third  onset  is 
also  difficult  to  quantify  because  of  the  spatial  separation  of 
the  stations.  On  the  basis  of  Figure  3  of  Tanskanen  el  al. 
[1987]  we  estimate  the  surge  to  be  between  Ivalo  ( L  =  5.7) 
and  Bjbrnbya  (L  =  9.4).  which  gives  a  distance  of  689  km. 


This  implies  a  Pi  2  period  of  312  s.  which  is  still  comparable 
to  the  measured  value  of  roughly  190  +  30  s  seen  12  min  after 
onset  in  Figure  9  The  important  thing  to  note  from  Figure  9 
is  the  definite  lengthening  of  the  Pi  2  period  as  the  surge 
expands  poleward  This  elTect  is  predicted  by  the  present 
model,  as  the  period  of  the  fastest  growing  mode  is  pro¬ 
portional  to  /.x. 

The  model  can  be  used  to  extract  other  important  infor¬ 
mation  about  the  onsets  using  this  important  data  set.  From 
the  value  or  \xh  of  4  4  km  s  and  noting  that  the  measured 
value  of  E„  was  approximately  20  mV  m  during  the  first  two 
onsets,  we  use  (3<i)  to  estimate  the  precipitating  electron 
energy  that  is  responsible  for  the  surge  movement  as  being  1.7 
keV.  (Note  that  Qll  =  5c1  2  ions  per  incident  electron  (see 
paper  1).  where  r.  is  the  precipitating  energy.)  Using  the  direc¬ 
tion  of  the  surge  motion,  one  can  also  estimate  the  magnitude 
of  the  north-south  polarization  electric  field  inside  the  surge. 
(The  measured  value  of  the  north-south  electric  field  compo¬ 
nent  is  25  mV  m.  which  includes  the  convection  electric  field 
which  may  dominate.l  Equation  (12)  of  paper  I  can  be  rewrit¬ 
ten  as 


tan  y,  -  (R  -  Er'En),H  1  +  REpiE0) 


(7) 
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where  ys  is  the  angle  of  the  surge  motion  relative  to  west.  R  is 
the  ratio  of  the  height-integrated  Hall  to  Pedersen  conduc¬ 
tivities.  and  (which  is  positive,  pointing  south)  is  the  polar¬ 
ization  electric  field.  (The  relation  Ep  =  £„/?(  1  -  a)  has  been 
used  in  deriving  (7).)  Equation  (7)  can  be  reexpressed  as 

Ep  E0  =  (R  —  tan  •/,)  (!  +  R  tan  /,)  (8) 

Grafe  [1982]  has  noted  that  R  ordinarily  ranges  between  1 
and  4.  with  some  rare  cases  in  which  R  is  greater  than  4.  It  is 
also  noted  that  R  tends  to  decrease  with  enhanced  substorm 
activity  and  widening  of  the  electrojet.  Using  (8),  it  is  found 
that  for  the  first  onset.  —2.2  mV  m  (R  =  I)  <  Ep  <  +2.1 
mV  m  (R  =  4).  which  is  consistent  with  full  closure  of  the 
poleward  Hall  current  into  the  magnetosphere.  (A  value  of 
2.75  was  used  for  tan  -/,.)  The  second  onset  is  like  the  first,  and 
the  third  is  too  complex  to  be  treated  here. 

We  conclude  therefore  that  the  three  onsets  discussed  are 
consistent  with  our  model.  Not  only  is  there  quantitative 
agreement  with  the  predicted  Pi  2  periods,  but  the  observed 
time  delay  for  the  Pi  2  pulsations  to  approach  a  monochro¬ 
matic  signal  is  compatible  with  a  frequency-dependent  decay 
time,  as  derived  above.  In  particular,  the  predicted  scaling  of 
the  Pi  2  periods  with  the  surge  size  and  poleward  velocities 
seems  convincing.  The  model  may  also  be  useful  as  a  tool  in 
separating  out  the  polarization  and  convection  electric  fields. 

6.  Discission 

We  have  presented  a  unified  model  for  the  motion  of  the 
westward  traveling  surge  and  the  generation  of  Pi  2  pulsations 
during  substorm  onsets.  Using  the  Inhester-Baumjohann 
ionospheric  current  system  shown  in  Figure  1.  it  was  shown 
that  the  zero-order  electron  precipitation  current  (del  con¬ 
trolled  the  surge  motion.  The  dynamics  of  the  surge  boundary 
was  found  for  a  constant  conductivity  profile,  which  led  to  a 
constant  value  of  the  boundary  propagation  speed  F.  This 
value  is  assumed  to  give  a  reasonable  estimate  of  the  initial 
surge  velocity,  which  in  a  more  realistic  treatment  must  be 
time  dependent,  as  the  conductivity  profile  evolves  in  time. 
The  first-order  equations  led  to  dispersion  equations  in  the  Pi 
2  frequency  range.  The  present  Pi  2  model  was  determined  to 
be  valid  between  surge  jumps  when  the  surge  was  considered 
quasi-stationary  or  when  the  surge  moved  sufficiently  slowly. 
During  jumps  it  is  suggested  that  a  nonlinear  theory  such  as 
that  of  Lysak  [1986]  is  more  realistic,  since  it  predicts  the 
presence  of  PiB  pulsations,  which  have  been  observed  by  Opye- 
noorlh  el  al.  [1980]  and  Bosinyer  el  al.  [1981] 

The  present  model  leads  in  a  natural  way  to  using  the  re¬ 
sults  of  Fridman  and  Lemaire  [1980]  and  Fdlthammar  [1978] 
in  determining  the  dc  precipitation  current,  which  plays  a 
major  role  in  damping  the  Pi  2  pulsations.  The  physical  re¬ 
quirement  that  the  Pi  2  pulsation  be  dampled  puts  a  lower 
bound  on  the  ratio  of  the  north-south  surge  dimension  and 
the  closure  parameter  x.  This  bound  depends  on  the  electron 
temperature  anisotropy  in  the  plasma  sheet  and  plays  a  cru¬ 
cial  role  in  determining  the  direction  of  surge  motion.  Grow¬ 
ing  models  were  found  to  be  quenched  by  the  poleward  ex¬ 
pansion  of  the  surge. 

Now  because  of  the  dispersive  nature  of  the  Pi  2  pulsation, 
its  frequency  composition  changes  with  time.  The  time  profile 
of  the  pulsation  also  is  very  sensitive  to  the  initial  conditions. 
If  all  the  modes  arc  excited  at  exactly  the  same  time  (i.e..  the 
same  phase),  then  the  time  profile  is  that  of  a  highly  damped 


wave  followed  by  periodic  sinusoidal  bursts  Even  a  small 
admixture  of  phases  al  I  =  0.  however,  will  cause  the  pulse  to 
exhibit  nonsinusoidal  behavior  just  after  onset  which  evolves 
into  sinusoidal  behavior  several  minutes  later  This  elVect  ap¬ 
pears  to  be  real,  as  seen  from  Figure  9  and  from  examples  of 
the  AFGL  magnetometer  data  (paper  21. 

Three  onsets  that  were  observed  on  June  22.  1979.  were 
compared  with  the  model.  It  was  found  that  the  poleward 
surge  velocities  and  the  estimated  north-south  surge  dimen¬ 
sions  led  to  predicted  Pi  2  pulsation  periods  that  were  consis¬ 
tent  w  ith  the  measured  v  dues.  The  period  of  the  observed  Pi  2 
pulsations  became  longer  as  the  surge  expanded,  which  is  con¬ 
sistent  with  the  model.  It  was  also  shown  how  the  model 
might  be  used  as  a  tool  in  determining  the  polarization  elec¬ 
tric  field  in  the  presence  of  a  strong  convection  electric  field 
Let  us  now  consider  a  scenario  for  auroral  breakup  based 
on  the  present  model.  Initially,  an  incipient  surge  region  of  a 
few  tens  of  kilometers  forms  near  local  midnight  [Akasofu. 
1974],  It  consists  of  an  enhanced  ionization  region  in  the 
ionosphere  that  contains  the  current  system  shown  in  Figure 
I.  Now  if  the  ionization  rate  at  the  surge's  poleward  boundary 
is  equal  to  the  ionization  rate  in  the  surge  interior,  the  surge  is 
stationary  [  Roiliwell  el  al..  1988],  If  the  rate  at  the  boundary 
is  greater,  then  the  surge  moves  poleward. 

It  is  interesting  to  compare  the  present  model  with  the  ex¬ 
perimental  observations  of  substorm  breakup  described  by 
Ikasol'n  [  1974.  p.  647], 

The  first  indication  of  a  substorm  is  a  sudden  brightening  of 
one  of  the  quiet  ares  (discrete  auroras)  lying  in  the  midnight 
sector  or  ihe  sudden  formation  of  an  arc  in  that  region.  A  typical 
subsiorm  develops  when  this  arc  is  located  near  the  equatorward 
boundary  of  the  belt  of  discrete  auroras,  but  near  the  poleward 
boundary  of  the  belt  of  diffuse  auroras.  In  most  cases,  the  bright¬ 
ening  of  an  arc  or  the  formation  of  an  arc  is  followed  by  its  rapid 
poleward  motion,  resulting  in  an  'auroral  bulge'  around  the  mid¬ 
night  sector  The  so-called  break-up'  occurs  in  the  bulge:  a  quiet 
curtain-like  form  appears  lo  be  disiupted  and  scattered  over  the 
sky 

In  terms  of  our  model  and  Akasofu's  observations,  auroral 
breakup  occurs  when  an  Inhester-Baumjohann  type  iono¬ 
spheric  current  system  forms  near  the  poleward  boundary  of 
the  diffuse  aurora.  The  poleward  Hall  current  in  this  embry¬ 
onic  system  is  closed  oil'  into  the  magnetosphere  through  the 
formation  of  the  arc  described  by  Akasofu. 

If  a  field-aligned  potential  drop  should  switch  on  over  the 
poleward  boundary,  the  arc  there  would  brighten  and  move 
poleward  with  a  speed  proportional  to  the  magnitude  of  the 
potential  drop  (see  papers  I  and  2  and  Rothwell  el  al.  [  1988]). 
The  boundary  motion  excites  standing  ionospheric  waves  in 
the  surge  interior  which  have  phase  velocities  which  arc  deter¬ 
mined  by  the  initial  boundary  speed  Because  of  their  long 
wavelengths  (  -  100  kn.)  these  ionospheric  waves  could  be  dif¬ 
ficult  to  detect  owing  to  the  turbulent  nature  of  the  ionosphere 
over  much  smaller  distances.  If  the  associated  frequencies  of 
the  ionospheric  waves  are  greater  than  or  equal  to  the  fre¬ 
quency  of  the  fundamental  toroidal  mode  of  the  attached  field 
line  [S/nc/er  el  al..  1981],  then  a  resonance  can  occur  through 
the  feedback  instability. 

The  resonance  is  considered  to  be  damped  by  electron-ion 
recombination  in  the  ionosphere.  Now  the  precipitating  flux, 
which  causes  the  damping,  is  a  function  of  the  electron  tem¬ 
perature  asymmetry  in  the  plasma  sheet  [ Fridman  and  Le¬ 
maire.  1980],  If  this  flux  is  insufficient  to  fully  close  off  the 
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poleward  Hall  current  in  the  ionosphere,  then  the  x  parameter 
is  less  than  1.  This  means  that  the  surge  will  move  predomi¬ 
nantly  westward  (paper  1)  and  the  coincident  Pi  2  pulsations 
will  be  strongly  damped.  (For  a  given  surge  size,  if  i  is  less 
than  1,  it  is  more  probable  that  the  threshold  criterion  for 
damping  is  met.)  A  value  of  a  less  than  1  implies  weaker 
precipitation,  which  for  a  given  field-aligned  potential  drop 
implies  enhanced  perpendicular  electron  heating  in  the  plasma 
sheet  according  to  the  model  of  Fridman  and  Lemaire  [1980], 
Perpendicular  electron  heating  in  the  plasma  sheet  would 
therefore  cause  predominantly  westward  surge  motion  and 
strongly  dampled  Pi  2  pulsations.  On  the  other  hand,  parallel 
electron  heating  in  the  plasma  sheet  (a  greater  than  or  equal 
to  I)  would  cause  predominantly  poleward  motion  and,  possi¬ 
bly.  the  momentary  presence  of  growing  Pi  2  pulsations.  As 
mentioned  above,  the  poleward  expansion  of  the  surge  limits 
the  amplitude  of  the  pulsations  when  the  damping  threshold  is 
reached. 

Therefore  according  to  the  present  model  the  substorm  is 
associated  with  enhanced  precipitation,  poleward  surge 
motion,  and  the  generation  of  finite  amplitude  Pi  2  pulsations. 
All  these  effects  are  common  features  associated  with  sub¬ 
storms  and,  as  shown,  may  be  caused  by  preferential  electron 
heating  in  the  plasma  sheet  parallel  to  the  magnetic  field. 

While  the  present  model  appears  to  be  very  successful  in 
explaining  a  variety  of  observations  during  substorm  onsets,  it 
is  by  no  means  all-inclusive.  It  is  still  a  linear  approximation 
to  a  highly  nonlinear  problem  where  the  nonlinear  effects  ap¬ 
parently  are  strongly  damped  over  relatively  short  time  scales. 
We  have  not  addressed,  for  example,  how  the  initial  Alfven 
wave  transmits  the  information  about  the  short  circuiting  of 
the  tail  current  between  the  magnetosphere  and  the  iono¬ 
sphere  during  the  creation  of  the  subslorm  current  wedge 
[Baumjohann  and  Glassmeier.  1984],  nor  have  we  touched 
upon  the  possibility  of  the  Alfven  waves  creating  a  field- 
aligned  potential  drop  through  turbulence  [ Lysak  and  Dum. 
1983;  Lysak.  1985].  Moreover,  a  more  realistic  treatment  of 
Alfven  wave  propagation  in  the  magnetosphere  and  ohmic 
losses  in  the  ionosphere  [ Goert :  and  Boswell.  1979]  is  needed. 
Finally,  there  may  be  more  than  one  source  of  Pi  2  pulsations. 
Edwin  et  aj  [1986],  for  example,  suggest  that  Pi  2  pulsations 
arise  from  an  impulsive  source  of  MHD  wave  energy  in  the 
plasma  sheet.  The  fast  magnetoacoustic  waves  that  are  pro¬ 
duced  propagate  dispersively  through  the  plasma  sheet,  which 
gives  rise  to  finite  Pi  2  type  wave  packets  nearer  the  Earth. 

It  is  clear  that  an  important  future  consideration  would 
involve  placing  the  mechanisms  developed  here  into  a  more 
global  framework.  In  this  regard  it  should  be  mentioned  that 
Zhu  and  Kan  [1987]  have  combined  the  temporal  propaga¬ 
tion  mechanisms  of  paper  I  with  the  global  assumption  of  the 
closure  parameter  developed  by  Kan  el  al.  [1984]  and  Kan 
and  Kamide  [1985],  Moreover.  Kan  and  Sun  [1985]  have  de¬ 
scribed  a  global  model  for  simulating  the  westward  surge  and 
Pi  2  pulsations.  In  this  model  they  combine  steady  state  as¬ 
sumptions  inherent  in  the  work  of  Kan  and  Kamide  [1985] 
with  a  description  of  Alfven  waves  bouncing  between  a  mag- 
netospheric  source  and  the  ionosphere.  Their  procedure  re¬ 
sults  in  a  temporal  scheme  involving  the  discrete  bouncing  of 
the  Alfven  waves  into  the  ionosphere  and  the  generation  of 
"steplike”  Pi  2  waveforms.  Note  that  in  both  works  the  iono¬ 
sphere  plays  an  active  role  and  is  a  source  of  Pi  2  pulsations. 
In  the  present  work  a  feedback  instability  is  the  driving  mech¬ 
anism.  while  in  the  work  of  Kan  and  Sun  [1985]  a  feedback 


instability  was  not  included.  Multipoint  measurements  are 
necessary  in  order  to  determine  which  of  these  models  are 
correct  under  dilferent  substorm  conditions. 


Am  NDIX 

Comparison  of  our  model  predictions  for  Pi  2  periods  with 
data  requires  knowledge  of  the  precipitation  energy  and  flux 
both  on  the  poleward  boundary  and  in  the  surge  interior.  This 
would  allow  a  direct  calculation  of  Vxi.  Since  this  information 
was  not  available,  it  was  necessary  to  assume  that  the  mea¬ 
sured  velocity  I'  of  the  poleward  boundary  was  equal  to  Ftl. 
Now  (3c)  implies  that  F  -  fV  only  when  the  ionization  rate 
along  the  poleward  boundary  is  much  greater  than  the  ioniza¬ 
tion  rate  in  the  surge  interior.  From  (3a)  we  see  that  Vxi  =  F<(, 
only  if  the  electron  precipitation  energy  in  the  surge  interior  is 
equal  to  the  precipitation  energy  on  the  surge  boundary.  We 
assume  that  this  latter  condition  is  satisfied  during  the  onsets 
considered  and  develop  a  model  to  ascertain  the  validity  of 
assuming  R,  «  Rh.  Now  enhanced  precipitation  on  the  pole- 
ward  boundary  implies  a  depiction  of  the  plasma  associated 
with  attached  llux  tubes  [Atkinson.  1984],  Therefore  as  the 
plasma  eonvects  earthward  through  the  surge  interior,  it  is  a 
weaker  source  of  precipitation  than  it  is  on  the  boundary,  and 
l  ir  F„  from  (3c). 

We  now  make  a  very  simple  model  that  illustrates  the 
physical  concepts  involved.  It  will  be  shown  that  for  reason¬ 
able  substorm  values  of  precipitation  current  (~  10  fiA  m2) 
along  the  poleward  boundary  the  precipitation  flux  in  the 
surge  interior  will  be  sufficiently  small  that  F  =  Vxl  from  (3c). 
This  will  validate  our  assumptions.  Consider  the  equatorial 
plane  and  a  plasma  with  a  number  density  ,V  =  0.3  cm  1 
conceding  toward  the  Earth  at  a  speed  F'  ^  50  km  s  [Huany 
and  Frank.  1986].  This  plasma  is  estimated  to  contribute  to 
the  electron  precipitation  fiux  over  a  field  line  length  of  ap¬ 
proximately  I.Rf  3.  where  I.  is  the  L  shell  on  which  the  pre¬ 
cipitation  takes  place  and  R,  —  6.37  x  lO*  m.  Therefore  the 
total  particle  fiux  which  is  available  for  precipitation  is 
N psl'cLRt/.„.  3.  where  /.ye  is  the  east-west  extent  of  the  surge 
region  as  seen  in  the  equatorial  plane.  This  flux  is  assumed  to 
be  totally  depleted  by  precipitation  along  the  surge's  poleward 
boundary  and  in  the  surge  interior.  The  total  electron  fiux 
precipitated  along  the  boundary  is  given  by  llll,cr-hcrve  e  and  in 
the  interior  by  jttie/xty  xr  e.  Here  7.N.  is  the  extent  of  the  con¬ 
ductivity  gradient  along  the  poleward  surge  boundary  (as  seen 
at  the  equator),  and  /.X1.  is  the  north-south  extent  of  surge  (also 
as  seen  at  the  equator).  Now  particle  flux  conservation  along 
the  field  line  implies  juhr/hls.y,.  =  jnh/.h/.f  and  =  .I'n/V-,.. 

where  the  quantities  without  the  subscript  e  are  defined  at  the 
ionosphere. 

By  equating  the  total  inward  convecting  fiux  to  the  flux  at 
the  poleward  boundary  plus  that  in  the  surge  interior  and 
using  fiux  conservation  along  the  magnetic  field  lines  we  final¬ 
ly  have  an  expression  for  jHh: 


Jus  ~ 


eNrxVR^  2 


-'/s!(/ii, >iisKz,  /.„)  +  i; 


A  m2 


(Al) 


where  the  azimuthal  scaling  factor  for  a  dipole  field  (/.  /A  = 
L3  2  [Lolkn  el  al..  1987])  has  been  used.  Now  if  we  make  the 
reasonable  assumption  that  7.,,  =  20  km  and  that  L  =  6  and 
f.x  —  280  km.  which  corresponds  to  the  first  two  onsets  of  June 
23.  1979.  then  (Al)  can  be  rewritten  using  the  numerical  values 
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given  above  as 

./'in.  =  22.5  J  14/n,  7M„  +  1 J  /iA  m2  (A2) 

Now  from  (3r)  (he  surge  is  stationary  if  the  current  ratio  is 
equal  to  1.  This  corresponds  to  a  boundary  precipitation  cur¬ 
rent  density  of  1.5  //A  m2.  For  V  ^  Vxi  the  current  ratio  must 
be  approximately  0.1  or  less.  This  implies  from  (A2)  that  jl{h 
should  lie  between  9.4  and  22.5  /(A  m2,  which  is  easily  real¬ 
ized.  We  therefore  conclude,  on  the  basis  of  this  very  simple 
convection  model,  that  the  assumptions  used  to  analyze  the 
Scandinavian  data  were  reasonable  and  appropriate. 
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